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Abstract 

We study polarization from scattering of light on a cloud in radial motion along the symmetry axis 
of an accretion disk. Radiation drag from the disk and gravitational attraction of the central black hole 
are taken into account, as well as the effect of the cloud cooling in the radiation field. This provides 
us with a self-consistent toy-model for predicted lightcurves, including the linear polarization that arises 
from the scattering. Strong gravitational leasing creates indirect images; these are formed by photons that 
originate from the disk, get backscattered onto the photon circular orbit and eventually redirected towards 
an observer. Under suitable geometrical conditions the indirect photons may visibly influence the resulting 
magnitude of polarization and light-curve profiles. Relevant targets are black holes in active galactic nuclei 
and stellar-mass Galactic black-holes exhibiting episodic accretion/ejection events. 
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1. Introduction 

The present-day evidence for black holes relies almost 
entirely on information carried by electromagnetic waves. 
X-rays play a particular role (Seward & Charles 1995): 
they are supposed to emerge from gas near a black hole 
horizon and bring us imprints of physical processes and 
conditions in the place of their origin and along the ray 
path. Various spectral and lightcurve patterns have been 
identified as likely signatures of supermassive black holes 
in galactic nuclei and stellar-mass black holes. These fea¬ 
tures presumably arise when matter is accreted from an 
immediate vicinity of the black hole. The ultimate goal 
of this effort is to prove the existence of event horizons in 
‘real’ nature and this way to discover black holes. Clearly, 
an affirmative proof is a great challenge that may still be 
far ahead of us. Here we study a related task, which ap¬ 
pears to be somewhat easier on the technical side, though 
it also represents an unresolved issue as yet: searching for 
rays of photons that encircle an ultra-compact star or a 
black hole. 

Accretion disks represent a common way of feeding 
black holes and generating photons, which then experience 
strong-gravity while travelling to a distant observer (Kato 
et al. 1998; Krolik 1999). Outside the disk plane density is 
much less, however, even regions near axis are not empty: 
winds and jets emerge, roughly along the symmetry axis. 
We address a question whether future observations of po¬ 
larized time-dependent signal can help recognizing the ef¬ 
fects of strong-gravity light bending, and what features are 
expected in light curves. To this aim we examine a model 
of warm clouds moving radially along the symmetry axis. 
Primary photons from the disk are Thomson-scattered 
and polarized. We assume that the process takes place 
near a Schwarzschild black hole, where the higher-order 


(indirect) light rays contribute to the observed radiation 
flux. The resulting modulation of intensity and polar¬ 
ization magnitude can reach a non-negligible level under 
suitable geometrical alignment of the black hole, the cloud 
and the observer. We treat the interaction of the cloud 
with the radiation and gravitational fields in the relativis¬ 
tic framework (Abramowicz et al. 1990; Vokrouhlicky & 
Karas 1991; Keane et al. 2001). This approach provides us 
with a self-consistent description of the cloud motion and 
the resulting observed signal. Both ingredients are conve¬ 
niently expressed in terms of the radiation stress tensor. 

There is the evidence for jets being formed only a 
few tens gravitational radii from a supermassive central 
black hole (Junor et al. 1999). The emission mecha¬ 
nisms producing the observed high-energy photons (X- 
and 7 -rays) are likely non-thermal, but it is not clear 
whether the synchrotron emission or the inverse Compton 
emission dominates in each particular source (Harris & 
Krawczynski 2002). Here we concentrate on the latter 
case (within the Thomson approximation), which seems 
to be relevant for radio-quiet sources with the ambient 
radiation acting on particles and fluids and determining 
their terminal (equilibrium) speed; see Noerdlinger 1974; 
O’Dell 1981; Sikora & Wilson 1981; and Phinney 1982 for 
original papers. This dynamical influence of the radia¬ 
tion held has been studied further by Sikora et al. (1996), 
Renaud & Henri (1998), Tajima & Fukue (1998), Fukue 
(1999) and Watarai & Fukue (1999). Ghisellini et al. 
(2004) proposed a model of aborted (failed) jets in which 
colliding clouds and shells occur very near a black hole 
and are embedded in strong radiation. According to their 
scheme, most of energy dissipation should take place on 
the symmetry axis of an accretion disc where the indi¬ 
vidual clouds either move away from the black hole, or 
they fall back. The process of gravitational and radiative 
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acceleration of plasma was studied also by Fukue et al. 
(2001) and Fukue (2005), who examined the efficiency of 
collimation towards the disk axis and applied this model 
to the case of microquasars. 

The effects of strong gravity are significant in this region 
and the challenge for future techniques is to identify sub¬ 
tle, yet specific patterns in X-ray lightcurves. Polarization 
studies could help to achieve this goal. Distinct features 
should arise from multiple trajectories of light rays con¬ 
necting the source with the observer along several different 
paths. The rays winding up around the photon circular 
orbit should experience a characteristic mutual time delay. 
We discuss the expected features in this paper. 

Polarization of light from scattering in winds and jets 
was examined by various authors. Following the early pa¬ 
pers (e.g. Dolan 1967; Angel 1969; Bonometto et al. 1970), 
Begelman & Sikora (1987) studied the linear polarization 
of initially unpolarized soft radiation up-scattered by cold 
electrons in a jet. Beloborodov (1998) examined the case 
of fast winds outflowing from an accretion disk slab. He 
found that polarization direction depends on the wind ve¬ 
locity; the terminal speed of the outflow plays a critical 
role. Poutanen (1994) and Celotti & Matt (1994) consid¬ 
ered the synchrotron self-Compton mechanism, a likely 
process producing polarization wherever magnetic fields 
interact with relativistic particles. The effect of electron 
temperature was also discussed: it reduces the final mag¬ 
nitude of polarization. Recently, Lazzatti et al. (2004) fur¬ 
ther discussed the Compton drag as a conceivable mech¬ 
anism for polarization in gamma-ray bursts. Horak & 
Karas (2005) studied the polarization of scattered light 
from a compact star, taking into account the light-bending 
effect. In fact, it was demonstrated that retro-lensing im¬ 
ages, which clearly require strong gravity, can give rise to 
specific polarimetric signatures in predicted lightcurves. 
Here we develop this model further by considering non- 
negligible temperature of the scattering medium and by 
changing the geometry of the primary source. We assume 
a Keplerian disk as the source, so we are able to examine 
situations that are relevant for accreting black holes. 

2. The model setup 

As a toy-model for non-uniform outflows and aborted 
jets we consider a cloud of particles moves through the 
radiation field of a standard thin accretion disc (Shakura 
& Sunyaev 1973). The disk defines the equatorial plane 
0 = 7r/2 of the system. Primary photons from the disc 
are scattered by electrons in the cloud, they are beamed 
(preferentially in the direction of the scatterer motion), 
and polarized by Thomson mechanism. We adopt the 
first scattering approximation (small optical depth r <§; 1 
of the cloud is assumed) and we restrict the cloud motion 
to the axis of symmetry, z = z{t). The basic setup of the 
model is captured in figure 1. 

The spacetime metric is (Misner et al. 1973) 

ds^ = (dd^ -b sin^ 6»d<()2), (1) 

where {t,r,0,(/)} are Schwarzschild coordinates, ,f(r) = 1 — 
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Black Hole 


Fig. 1. Geometry of the model. An accretion disk is the 
source of primary unpolarized light, which is then Thomson 
scattered on a cloud. The cloud moves radially along the 
axis, 2: = z{t), in the radiation field of the disk. Gravity of 
the central black hole infiuences the motion of the cloud. The 
light rays of primary and scattered photons are also affected. 
Direct and indirect (retro-lensed) light rays have different de¬ 
gree of linear polarization and they experience different am¬ 
plification and the Doppler boosting. The observed lightcurve 
is produced by all the rays reaching an observer at view an¬ 
gle i far from the centre (along Z-direction). The indirect 
photons contribute most significantly to the the total signal 
if the cloud moves toward the black hole and the observer 
inclination is small. 

Rs/r is the redshift function, Rs = "IGMjt?. Geometrized 
units c = G = 1 will be used. 

The total bolometric luminosity of a standard disk, L = 
GMM/2i?in, can be scaled by the Eddington luminosity 
Le = dncGMimp/ We thus define A = L/Le- The 
corresponding intensity emitted from the disk is 



where we denoted, by an asterisk, the intensity rescaled 
to the dimension-less form. Hereafter the radiation 
terms will be made dimension-less by scaling them with 

{nipC^/ a^Rs) K. 

Photons travel along null geodesics and give rise to the 
radiation stress tensor at every point outside the disk 
plane and above the black hole horizon. Figure 2 shows 
the dependence of the components relevant for our cal- 
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Q = pYY_pXX^ ( 4 ) 

U = -2P^^, (5) 

where bars denote quantities in the comoving orthonormal 
polarization frame, X-axis of which is given by projecting 
the symmetry axis onto the observing plane and Z-axis co¬ 
inciding with the direction of scattered photons. are 

local components of the radiation stress tensor (the spatial 
part of the energy-momentum tensor E = Tbllb 

is local energy-density. The values of Stokes parameters 
(3)-(4) are normalized hy A = 3T/167r, r = rieaTR is the 
Thomson optical depth of the cloud in terms of electron 
density rie, i? <C i?s is radius of the cloud. Polarization 
is purely linear (V = 0) and its the magnitude is equal to 

n = (Q2 + i/2)i/2/j_ 


Fig. 2. Normalized components of the radiation stress ten¬ 
sor on the symmetry axis for radiation originating from a 
standard disk. The inner edge of the disk corresponds to the 
innermost stable circular orbit, = Jtisco = 3its; no in¬ 
tervening matter is considered below that radius. The outer 
edge is at Rout = 10®its- Location of the photon circular orbit 
i?ph = 1.5i?s is also indicated. 


culation on symmetry axis of the black hole - accretion 
disk system (we employ similar notation for the radiation 
stress tensor components as e.g. Fukue & Hachiya 1999; 
Watarai & Fukue 1999). 

In spite of the fact that we employ a particular model of 
a thin accretion disk in this paper, several rather general 
features are captured in Fig. 2 that we expect to shape 
the radiation field of astrophysically realistic models too. 
These can be summarized as follows. 

— The radiation stress tensor vanishes for z —> oo (^ —> 1) 

because the angular size of the disk, as observed by 
distant observers, decreases to zero. 

— Far from the centre, te-component of the stress tensor 

(measured by a static observer) is positive, however, 
this term changes its sign at certain value of height z 
near above the disk. This effect operates for a static 
observer inside photon circular orbit, who intercepts 
photons falling into the black hole. 

— All three components diverge as z ^ Rs (^ ^ 0). Close 

to the horizon the accretion disk appears again as 
a point-like source producing the energy density oc 
^“^(z). Here, two effects act against each other: 
the gravitational bending of light-rays reduces the 
solid angle that the disk occupies on the observer 
local sky by factor oc p whereas the gravitational 
redshift increases the radiation intensity by factor 
a (blueshift occurs, in fact). 

Intensity of scattered light and its polarization are 
characterized by Stokes parameters (Chandrasekhar 1960; 
Rybicky & Lightman 1979). These can be obtained by 
integrating over all directions n; of photons incident on 
scattering particles of the cloud from different points in 
the disk (Horak & Karas 2005): 

I = E + P^^, 


3. Polarization from scattering on warm clouds 
near a black hole 


So far we assumed the cloud to be cold and we did not 
consider any microscopic (random) motion of the scat¬ 
tering electrons within the cloud volume. Hereafter, in 
order to account for the electron temperature we aver¬ 
age the equations (3)-(5) over the electron distribution 
n{l3g) = Uefiye) in the cloud comoving frame (CCF), 
where the electron distribution is entirely isotropic and 
7 e = (1 — /3e)~^^^ is the corresponding Lorentz factor. 
Stokes parameters are first evaluated in the local electron 
comoving frame (ECF), Lorentz transformed to the CCF, 
and the result is averaged over random velocities of the 
electrons. We obtain 


I={1 + A){E + P^^) +B{E- 3P^^) - 2AP^, (6) 

Q = pYY_pXX^ (7) 

U = -2P^^, ( 8 ) 

where P^ is the radiation flux component measured in 
direction of the photon after scattering. We introduced 
mean values 


A = 


— 4 

3 




(9) 


(x) = J X f {je) die denotes averaging over the electron 
energy distribution; see Horak (2005) for details of the 
derivation of eqs. ( 6 )-( 8 ). 

The Lorentz boost from CCF to the static frame intro¬ 
duces the familiar factor where 'D = i~^{l — /3cos'd)~^ 
is the Doppler factor, 7 = (1 —/3 is the cloud bulk 
velocity and is the angle between this velocity and the 
scattered photon direction: 

I = V^P Q = V'^Q, U = V'^U. (10) 


In this way the frequency-integrated Stokes parame¬ 
ters are found in terms of the radiation stress tensor 
components. The relation between the CCF compo¬ 
nents {E, P^^) and the static frame components 

{E,F‘^,P‘^^) is 


pXX _ pp. 


P^^ (cosi9 - (3f + y^Asin^ ^ 


( 11 ) 


( 3 ) 
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Fig. 3. Left: critical velocities of the cloud motion, /5i(^) and /32(0? which the observed polarization vector changes its orientation 
between the longitudinal (inside the shaded region) and the transversal one. The saturation curves /3 = /5o($) are also plotted and 
labelled by the values of the normalized luminosity A (for A = Acrit the saturation velocity vanishes at infinity). Right: several 
examples of possible trajectories of a warm electron cloud (thick curves). The initial condition is /3(to) = ^5 the values of ^(to) a.re 
indicated by a circle. A monoenergetic distribution of scattering electrons was assumed, /(7e) = ^(7e — (7e))- Cold cloud trajectories 
are shown for the sake of comparison (thin solid curves). The saturation curve is also plotted (dashed), corresponding to the 
luminosity parameter A = lOfi. 


pZZ _ jy2 


P^^sm^d + 'y^A{cos'd-pf , 


( 12 ) 


pYY^pyy^ (13) 

E ="f‘^{E-2(3F^ + (14) 

=Vj^{cos9-P)W. (15) 

where we denoted A = P'^E — 2/3F^ + P^^ and W = (1 + 
P^)F^ — P{E + P^^). Axial symmetry ensures that P^y 
vanishes for the primary radiation field originating from 
the disk and impinging on the cloud. 

Properties of the scattered signal in the static frame 
of course depend on the motion of the scatterer, which 
itself results from the interplay of acceleration mechanisms 
acting on the cloud. In our model the total four-force /“ 
is a superposition of the radiation and inertial terms. The 
radiation term is 


/“ = -a^UeSV \cT°‘Pui3 - (A + C)r^'"itpu„u“l, (16) 


where (5V is the cloud volume, C = 1 + ^A. The change 
of the cloud internal energy by cooling is 


^ (e^ - 2PF: + P^P:^) , (17) 

where we expressed proper time s in units of Rs/c and 
we denote /i = ~ 1 (for electron-proton plasma) or 

/i ~ 10^ (for electron-positron plasma). The parameterisa- 
tion by /r allows us to consider different types of plasma. 
However, we note that the radiation field of a standard 
disk alone is unable to efficiently accelerate protons, and 
so the case of electron-positron plasma seems to be a more 
relevant application here (Fukue 2005). 

Dynamics of the cloud as a whole is governed by equa¬ 
tion of motion in the form 


ds 


(7e)-'CAM {1 + P^)f: 


-/3(A* + Py) 


Rs 

27^71/2 ’ 


(18) 


which includes also the black hole gravity. The last two 
equations couple the cloud trajectory with the evolution 
of its internal temperature. In order to close the set of 
equations we need to assume a specific form of the electron 
energy distribution /(7e). 

In the case of cold clouds (A = 0, C = 1) the dynam¬ 
ics is characterized by the saturation velocity PoiO to 
which particles are asymptotically accelerated if inertial 
effects are small compared to the radiation force (Sikora 
& Wilson 1981; Icke 1989). This condition corresponds to 
the limit of A —> 00 , d/3/ds = 0. We thus obtain 


P{z) ^ Po = a- \/cr2 - 1, 


O' 

2Fp 


(19) 


Further, for polarization there are two critical velocities, 
/3i(^) and /32(C) i which the polarization vector changes 
its orientation between transversal and longitudinal one 
(this can be seen from the condition (5 = 0). Similar ef¬ 
fect of polarization direction changing with velocity of the 
scattering medium was studied by Beloborodov (1998; for 
scattering in fast outflows from accretion disks) and Horak 
& Karas (2005; for radially moving clouds near a compact 
star). This behaviour is demonstrated in the left panel 
of figure 3. The right panel shows the corresponding in¬ 
fluence of the radiation cooling of the cloud on its bulk 
motion. 

From eq. (18) we find that the motion of warm clouds 
is governed by the same equation as cold ones, provided 
that the normalized luminosity is rescaled: 

A^A' = {7e)-iCA. (20) 


In other words, by keeping the disk luminosity at a con¬ 
stant value of the Eddington parameter, the cooling (17) 
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Fig. 4. Comparison between two typical cases with the identical initial conditions except for the cloud temperature: a cold cloud 
(left panels, (ye) = 1) versus a warm cloud (right panels, (ye) = 3 at start). Examples are shown of the time behaviour of the observed 
intensity (top row) and the magnitude of polarization (middle row). Contributions of the retro-lensing images have been summed 
together and plotted (by a dashed line); they are clearly distinguished from the signal produced by the direct-image photons (solid 
line). The retro-lensing photons arrive with a certain time delay with respect to the direct photons. The delay is caused by photons 
taking different path and circling at r Cr: i?ph (the effect can be recognized by comparing the relative time shift of the features). We 
also show, in the bottom row, the corresponding velocity profile /3(^) of the cloud motion. The temperature evolution is given by 
solving eq. (17) along the trajectory: the warm cloud cools down to (ye) = 1.3 along its entire track. Polarization vanishes at the 
moment when the cloud crosses one of the curves /3i (^), 02 (^) • In both cases the view angle was 2 = 5 deg, the disk luminosity A = 2. 
The initial distance z{to) = (i.e. ^ = ^), /3(£o) = 0.8. 
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gives rise to trajectories gradually deviating from the cor¬ 
responding cold-cloud limit. 

The effect of changing temperature is such that the sat¬ 
uration velocity curve gradually moves across (/3,^)-plane. 
This helps us to understand the form of trajectories, which 
create loops (see fig. 3, right panel). The critical point is 
given by condition /3 = 0, 0<^<oo and it defines the 
distance where equilibrium between radiation and gravity 
can be reached. This is a saddle type point, character¬ 
istic to an unstable equilibrium, which makes the system 
behaviour qualitatively different from the case with a cen¬ 
tral star as a source of primary irradiation (cp. with the 
critical points examined in Abramowicz et al. 1990). Our 
result here resembles the case of clouds with non-constant 
mass, discussed by Keane et al. (2001). 

Now we can find the Stokes parameters of scattered light 
along the cloud path. When determining the temporal 
evolution of observed intensity and polarization we con¬ 
sider first three images in the observed radiation - the di¬ 
rect (zeroth-order) image and two first-order images. The 
latter are retro-lensing images formed by rays making a 
single round about the black hole by the angle = 27r±j, 
where i is observer inclination (we denote the retro-lensing 
images by ± signs to indicate that they wind up in the 
opposite directions around the black hole). The higher- 
order images give a progressively weaker contribution to 
the final signal, hence we can safely neglect them. The 
retro-lensing photons give rise to peaks in the observed 
signal occurring with a characteristic mutual time lag af¬ 
ter the direct-image photons. Duration of these features 
is very short - comparable to the light crossing time - 
and the time span of the plots can be therefore scaled by 
tic ^ 1-5 X 10-'‘M/(10 Mq) [sec]. 

Typical profiles of the resulting lightcurves and polar¬ 
ization curves are shown in figure 4. Here we notice the 
difference in the lightcurve profiles and the the polariza¬ 
tion curves, which is caused by effect of the cloud initial 
temperature and its gradual cooling. The contribution of 
the retro-lensing images to the total observed flux is max¬ 
imum when the scatterer, the black hole and the observer 
are well-aligned (i ^ 0) and the scatterer moves towards 
black hole (/3 —> — 1). In this situation a non-negligible 
fraction of photons are scattered toward on the light cir¬ 
cular photon orbit, their energy is Doppler boosted, and 
eventually they are redirected toward the observer. On 
the other hand, polarization magnitude attains maximum 
in a slightly different direction whereas, for symmetry rea¬ 
sons, n vanishes for a strictly aligned observer. The direct 
image photons are of course most intense in the opposite 
case (/3—> 1), and so there is an interplay of different orders 
of the images sensitive to the model parameters. 

Identifying the short-duration features in lightcurves 
could help confirming the existence of the photon circu¬ 
lar orbit around the central body, for which black holes 
are the most likely cause. Polarimetric resolution provides 
additional information about the source and it could be 
used to constrain the parameters of the black hole. 


4. Discussion and conclusions 

In X-rays, future polarimeters could be employed to 
probe jets and winds in strong gravitational fields of the 
central compact object. Polarimetry is a powerful tool 
that can provide additional information, which would be 
difficult to obtain by other techniques such as traditional 
photometry and spectroscopy. In this way polarimetry 
helps to discriminate between different geometries and 
physical states of sources where accretion processes are ac¬ 
companied by fast radial motion of the blobs of material. 
Naturally, this goal would require sufficient sensitivity in 
X-rays; the scattered signal is mixed with primary pho¬ 
tons, which reduce the final polarization. We have seen 
that the predicted features are flashing for only a brief 
period of time and the maximum polarization degree is 
typically ~ 1 percent or less. 

In order to allow the analytical treatment we employed 
various simplifications: we considered the bolometric 
quantities and the Thomson cross-section for scattering of 
primary photons on electrons (rather than Compton scat¬ 
tering and the Klein-Nishina cross-section; see Melia & 
Konigl 1989; Skibo et al. 1994; Madau & Thompson 2000). 
We also assumed that the flow is not magnetically dom¬ 
inated, although astrophysically realistic models require 
magnetohydrodynamic effects to be taken into account 
(Begelman et al. 1984; Beskin et al. 2004). Effects of 
general relativity were taken into account in the limit 
of the non-rotating black hole spacetime (we neglected 
the effects of frame-dragging for the sake of simplicity; 
cf. Vokrouhlicky & Karas 1991). Likewise, we adopted 
the simplest possible parameterization of the disk emis- 
sivity via the standard Shakura-Sunyaev model; this 
could be improved by including relativistic effects (Page 
& Thorne 1974) and a more realistic description of the 
disc itself. These changes will be necessary to provide 
quantitative and astrophysically realistic results for the 
polarization degree, however, we do not expect any quali¬ 
tative change regarding the signature of indirect photons. 
The predicted polarization is either parallel or perpendic¬ 
ular to the projection of cloud velocity onto the observing 
plane. 
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